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ABSTRACT 

This thesis attempts to develop a program for the Whirlwind 
computer operating in the interpretive mode for the flux distribution 
in a finite cylindrical reactoro 

The ieee chosen for investigation was a hypothetical organic 
moderated, highly enriched, cylindrical reactor having a power output 
capable of fulfilling one half the propulsive requirements of a 
nuclear powered naval cruisero 

Heat transfer and basic physics calculations were accomplished 
to determine preliminary size and material requirement of such a 
reac toro 

The flux equations were derived using numerical integration 


methods suitable for computer programmings 





Conclusion 


1) The code as presently programmed is incorrecte 

2) The grid spacing as used is much too coarse. This 
difficulty is not inherent in the coded programe 

3) The mode of computer operation is unsuitable for an 
iterative process such as this because it results in 


. a code that progresses too slowlyo 


Recommendations 

1) Prior to the use of any computer a formal or semi-formal 
course in programming should be completed. 

2) A grid spacing comparable to the mean free path of the 
neutron should be selected. 

3) The fastest possible mode of calculation on the computer 
available should be usedo 

l.) For a coded program of this difficulty a faster machine 


with a much larger fast memory is very desirable. 


Thesis Supervisor Melville Clark, Jr. 


Title Assistant Professor of Nuclear Ingineering 
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INTRODUCTION 


The desirability of a uniform power distribution (P wee =] 
has long been realized since it will not only decrease the mass of 
the fuel required for the attainment of criticality but also will 
increase the total power output by allowing the whole core instead 
of only a small fraction of the reactor to operate at the highest allow- 
able temperaturee Operation at the highest possible temperature allowed 
by metallurgical considerations over the entire reactor will increase 
the average core life since more uniform burn up is attained. 

It is possible to achieve this uniform power distribution by 
various meanSe One of the first methods to be used was that of 
poisons. In this method long lasting poisons are placed in the 
reactor structure in the anticipated locations of high flux. By 
placing poison in these regions the peak flux is reduced to such a 
level that the power density is approximately uniforme This pro- 
Geiure2is a very wasteful one from the point of view of neutron economy. 
Toe regions of the reactor that are contributing greatly to the overall 
neutron efficiency of the reactor are penalized to the extent that 
they can be no more efficient than the most inefficient region. For 
power producing reactors, particularly mobile ones, this waste of 
neutrons is very undesirable. 

The use of poison is the first of two methods that exp]lsit a 
nonuniform macroscopic cross section to achieve the uniform 


distribution of power. The second consists in the use of a nonuniform 


macroscopic fissigqn cross sections This method allows each region of 


a 





the reactor to operate most efficiently. The efficiency is determined 
by the position of the region within the reactor. The leakage from 

a reactor causes a reduction in the value of the flux at large 
distances from the center of the core. To compensate for this reduc- 
tion in flux the macroscopic fission cross section is increased in 

the regions far from the center and decreased in the central portions. 
To this end, the fuel concentration is increased near the edgese 

Neutron economy can also be improved by surrounding the core with 
neutron reflector material. A light substance having a relatively 
low absorption cross section performs this function very well, ic@e, 
heavy water or berylliume 

Since the reactor to be investigated is intended for a hypothetical 
naval cruiser, the use of poisons to achieve a uniform power distribu- 
tion is not desirable, for a compact reactor is requirede Though the 
reactor must be compact, it is desirable to install a reasonably thick 
reflector for the purpose of improving the neutron economy, attenuating 
“the-neutrons Leaking out, and cooling of the reactor shell and 
pressure vessel. 

It would be possible to have a nonuniform macroscopic cross 
section of the reactor in each of its directions. For the spherical 
reactor the variation would be radial only, but the application of a 
spherical reactor to shipboard use is impracticale The noriyal con- 
figuration of spaces aboard a ship is rectangular and a spherical 
reactor does not adapt very well to this arrangement. However, a 
rectanguiar reactor allows efficient utilization of space aboard a 


ship, but from the miclear point of view this shape is not too desir=- 
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able because leakage depends on the relationship of surface area to 
volume. It can be seen that the neutron economy decreases in 
progression from the sphere to the cylinder to the rectangular 
parallepipeds 

With these considerations in mind a right cylindrical reactor 
shape will be investigated. With this shape it would be possible 
to vary the properties in two directions, radially and axiallye 
While this goal is theoretically interesting, it is prohibited 
when cost is considered. The reactor that'is tobe investigated will have a 
uniform composition in the axial direction with finite radial segments 
differing in properties. 

For the shipboard use of a nuclear reactor power plant it is 
desirable to have a plant as light and safe as can possibly be 
deSigned within the space allowed but still without incorporating 
undue quantities of the "exotic" and expensive materials or coolants. 
One possible solution is the use of a highly enriched reactor using 
a coolant that permits high operating temperatures at low pressures 
and does not become dangerously radioactive while passing thru 
the reactor coréo 

A coolant that has a low vapor pressure at high operating 
temperature will allow a system pressure which is much lower than that 
employed in the present day water cooled and moderated shipboard 
reactor installation. The use of such a coolant will not only result 
in much better steam conditions, but it will also make possible an 
overall weight reduction, since a low pressure system can be usede 


The saving in weight is particularly significant in the pressure 
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vesseélo 

If the coolant does not become radioactive in,its passage thru 
the reactor core, the shielding can be concentrated primarily around 
the reactor proper, provided that impurities in the coolant are 
kept at a very Low concentration. Even with a small amount of 
impurities, the radiation emitted by the coolant would be low 
enough to permit reduced thicknesses of shielding around the primary 
LOOPe 

A liquid organic hydrocarbon will fulfill the requirements listed 
above, but also will present some unsolved problems. It will (1) tend 
to increase the size of the reactor by decreasing the hydrogen moderate 
tion relative to a water cooled and moderated reactor, (2) tend to 
increase the reaetor size due to its less efficient heat transfer 
capabilities compared with water and (3) exhibit both thermal and 
radiolytic breakdown in which the original hydrocarbon will, primarily 
by polymerization, change into longer chain hydrocarbons with the 
evolution of relatively large amounts of gas and a great increase in 
viscosity of the coolarte 

The latter problems will be investigated at the National Reactor 
Testing Station on the Organic-Moderated Reactor Experiment (ORME). 
This experiment is not a prototype for an organic power reactor but 
is being built to test the feasibility of its use as a reactor a 
At the West Milton Annex of the Knolls Atomic Power Laboratory the 
Atomie Energy Commission plans to investigate the feasibility of an 
organic moderated and cooled reactor for naval propulsione This will 


be conducted on the Naval Organic Reactor Experiment, (NORE )o 


4 





The heat transfer properties of the various liquid organic 
hydrocarbons that are presently under consideration are not generally 
available. Some presently under investigation are diphenyl, 
O-terphenyl, M-terphenyl and P-terphenyl. A proprietary heat transfer 
medium "Dowtherm" is available commercially and is composed of 7305% 
diphenyl oxide and 26.5% diphenyle It has a working range of }50- 
750°F. at pressures less than 15 psige ‘Since "“Dowtherm" closely 
resembles the hydrocarbons that will be investigated on the ORME and 
its heat transfer properties are readily available 14% will be used 
in this investigation. 

The fuel element chosen for this investigation is of the MIR 


[3] 


(Materials Testing Reactor) type. An effective length of element 


{his length was chosen so that the overall height 


was set at },6", 
of the reactor, control devices and shield would fit into the present 
day engineering spaces of a cruisers Any additional height would have 
an adverse effect on stability, since it would result in a higher 


shield deck and thus an increase in the height of the center of 


gravity of the shipo 





Fuel Element Dimensions 


t # 0.020" 
t = 0.015" 
Ww = 208 

L = 18.0" 


t = Q.115" (Case 1) 


= 0.0919" = 0.8 t (Case 2) 


Ss 
a 
s 


t 


t = 0,069" = 0.6 +t (Case 3) 
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PROCEDURE 
Heat Transfer ' 

The controlling variables in this study were (1) the maximum 
surface temperature of the fuel element or the maximum temperature 
at the centerline of the fuel element and (2) the maximum coolant 
speed. The maximum surface temperature of the fuel element may be 
controlling, since the preliminary investigations on the ORME have. . 
shown that above 800° F the thermal breakdown is increased eattag. 
However, liquid loading requirements of a naval vessel make storage 
of used coolant a simple matter and in the interest of higher 
efficiency some breakdown will be accepted. Therefore, a fuel element 
surface temperature less than 850° F will be accepted. If the maximum 
fuel temperature at the centerline is allowed to exceed a value of about 
1130° F, a phase transformation of uranium is encountered in which an 
anisotropic expansion and possible rupture of the fuel element may 
result. Therefore, a second design criterion will be required that 
the fuel element centerline temperature be less than 1130° F. 

Preliminary studies on the ORME have also shown that coolant 
speeds greater than 15 feet per second result in excessive erosione 
As a third design criterion, it will be required that coolant speed 
be less than 15 feet per seconde 

The first investigation was with a reactor coolant outlet temper= 
ature of 750° F with various temperature rises across the reactore 
This investigation was amplified by varying coolant speed while holding 
the temperature rise across the reactor constant. Each of these 
Situations were further studied by considering the MTR fuel element 
modified to the extent that the plate spacing was narrowed to 80% 
and 60%. 





' Nomenclature for Heat Transfer on pages 22 and 23. 
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Though it was known in advance that the speed would be kept at or 
less than 15 feet per second calculations were pursued with higher 
values. This was done to show that a reduction in size can be accom- 
plished in the futures if the erosion problem at high coolant speeds 
can be overcomes 

A heat rate requirement of 6.06 x 10° Btu/hr was chosen based on 
the full power steaming rate of a typical cruiser. 

With this value of Q (6.06 x 10°) a coolant outlet temperature, 
coolant temperature rise, coolant speed and plate spacing are chosenoe 
Since the inlet and outlet temperatures are known the average values 


of density and specific heat may be calculated. Formulae l, 2, 3,5 and 


are used respectively to calculate w, Ag G and no 











Since the dimensions of the ‘fuel element for each case are known, 
both the hydraulic diameter and the total flow area can be calculated. 


With the product of the mass flow rate and the hydraulic diameter the 


produce of the heat transfer coefficient and the hydraulic diameter can 


be read from the graph in the "Dowtherm" booklet. Formulae 5, 6, 7, 


¢ 


(1) 


(2) 


(3) 


(4) 





d 8 h i A AT 
an are then used to calculate, respectively, ae T fuel? “ead? 
and 06 


ee (5) 


a ek eee (6) 


=< lal ai (7) 


Q=<f -tT -—— (8) 


Datum from the calculation of a coolant outlet temperature of 750° F 
for various values of temperature rise and speed for cases 1, 2 and 3 
are in Tables I, II, and III and Figures I, II, and III respectively. 

The cladding temperature existing for each of these conditions is 
platted on Figure IV. 

The narrowest plate spacing, case 3, gives a significantly smaller 
reactor. Since the pumping power of any of these cases proved to be 
wery small compared with the total power output, it was decided to 
investigate the result of a decrease of coolant outlet temperature for 
this case only. Coolant outlet temperature was reduced to 700° V and the 
speed was held at 15 feet per second. Datum from this calculation are 
in Table IV. These data were plotted as a dotted line on Figure lVe 

For these twelve conditions the fraction of uranium in the uranium 
zirconium dispersed fuel plate was calculated and plotted on Figure V. 

It is desirable that the fuel be of such a composition that it contributes 
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TABLE I 


Datum for coolant outlet of 750°F 
and wide plate spacing 


A. REACTOR AT = 200°F 


V 


(ft/sec) (1b/sec-ft* ) 


5 
15 
25 
50 


a 
(ft/sec) (Btu 


15 


G 


23601 
709.0 
1180.0 
2360.0 


x 10 


1.86 


Be REACTOR AT = 150°F 


Vv 


(ft/sec) (1b/see-£t*) 


5 
15 
25 
50 


Vv 


1 


G 


23109 
70520 
1160,0 
2318.0 


Pes 108 
(ft/sec) ¢(Btu/nr-ft 


3.60 


8 


/nv-ft?) 


+ 
excessive temperature 


Yoranium content too low 


(case 1) 


h > Ap 
(Btu/hr-ft -°F) (ft?) (no. of plates) 
686. 5.36 
1611. 1-79 
2310. 1.07 
1,090. 0.5h 
rT, Tali, ol 
(Pr) (°F) (oF) (FF) 
2007 6262 215.8 1015.8* 


n 


2685 
896 
538 
269 


h ° Ap n 
(Btu/hnr-ft -°F) (£t°) (noe of 
682, 7026 35639 
1591. 202 15211 
2320. 1eh5 727 
1,00. 0.73 36), 
To-Ty Ty-To To-T, To cm 
(°F) "Cr) “@r) (or) % 
15eh 4600 19763 94703” 1965 


plates ) 





Table I (Contimed) 


Ce REACTOR AT = 100°F 


Vv G . h 9 Ap n 
(ft/sec) (1b/sec-ft) (Btu/nr-ft°-°F) (ft2) (noe of plates) 
i) 227.0 669. 11.01 55520 

15 681. 1569. 3.67 1,839 
25 1135. 2295. 2.02 1,103 
50 2270 1,000, 1.10 552 
v Q x 10°, ee oe a fa 
(ft/sec) (Btu/hr-ft') (°F) (°F) (°F) (°F) % 


15 2.37 10.1 30.3 133.8 883.8" 1395 


De. REACTOR AT = 50°F 


Vv G 5 h 9 Ap al 
(ft/sec) (1b/sec=-ft ) (Btu/hr-ft -°F) (£t2) (noe of plates) 
5 221.8 659. 22.10 11,230 
15 667.0 1556 Tel5 35730 
25 1410,.0 22755 he 7 2,210 
50 2218.0 3910. 2.2h 1,123 
8 
9 Q x 10 T.-T T-T. Teor T f 
v : li? tae 2 
(ft/sec) (Btu/nr-tt?) (er)’ (*F) (°F) (¢F) a 
15 Dea, he97 15-0 66.2 816.2 6099 


+ 
excessive temperature 


+ 
uranium content too low 


LO 





TABLE Il 


Datum for coolant outlet of 750°F, 
coolant speed of 15 ft/sec, and medium plate spacing 


+ 
temperature excessive 


uranium content too low 


(case 2) 
AT i, h , n Q_ x 10°, 
(°F) (lb/sec-ft ) (Btu/hr-ft -°F) (no. of plates) (Btu/hr-ft" ) 
200. 709. 1671. Tea 3.89 
150. 7056 1660. 15156 2.88 
100. 681. 16396 2300. 1.89 
50. 667 1631. 1660. 0.9) 
AT T.—0 T=T ‘Tee T f 
2 ob 2 u 

(°F) (*P) (°F) (°F ) (F) % 
200 1606 19.7 19309 9369 19.6 
150 mbes 36.9 Lh bro 2 89.0 16.1 
100 8.1 23 96.3 8603 11.0 

50 1.0 12,0 706 79706 eae 





TABLE IIL 


Datum for coolant outlet of 750°F, 


coolant speed of 15 ft/sec, and narrow plate spacing 


(case 3) 
AT G : h ; n x 10°, 
(°F) (lb/sec-ft ) (Btu/hr-ft -°F) (no. of plates) (Btu/hr-ft  ) 
200 709 1796 11:96 209 
150 705 1775 2020 2.2 
100 681 1731 3061 lel 
50 667 1718 6230 0.70 
AT { <T Tee T = T f 
oa 1 2 2 2 u 

(°F) () (°F) (°F) ( °F) % 
200 12h 3763 13502 885.2" 1603 
150 962 2706 101.3 851.63 12.) 
100 sha 18.2 68.5 818.5 Boh 
50 3.0 900 3309 78309 hod! 


*% 
temperature excessive 


y 
uranium content too low 
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TABLE IV.. 


Datum for coolant outlet of 700°F, 
coolant speed of 15 ft/sec, and narrow plate spacing 


(case 3) 

AT Pr. h , n Q_ x 10° 
(°F) (1b/sec-ft ) (Btu/hr-ft -°F) (noo of plates) (Btu/hr-ft" ) 
200 73905 1826 871 5.0 
150 726.8 1812 116), 307 
100 71265 1775 1767 205 

AT tt, de a fe 
(°F) (°F) (°F) (°F) (°F) % 

200 2103 6h. 0 228.) 92hel; 26.0 
150 15.9 h79 ieee 87262 20.3 
100 10.5 31.5 1556 815.8 1.0 


x 
temperature excessive 


‘6 
uranium content too low 
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to the structural strength of the reactor. A uranium content of 10% by 
weight or greater has desirable strength characteristics. 

Of the six conditions that have a fuel element surface temperature 
of less than 850° F only two have a uranium fraction of above 10%. 


These are shown in Table V, 


t T case n T 
co 2] 
750 100 5 3061 818 
700 100 3 1767 815 


Table V Particularly Interesting Cases 


Figure VI shows the variation in number of plates for case 3 
holding the coolant speed at 15 feet per second and varying the 
temperature rise across the reactor. 

The size resulting from the condition with a coolant outlet 
temperature of 700° F is better for the following reasonsse 

1) it is by far the smaller, 
2) the increase of uranium fraction by 2 1/2% gives a 
significant increase in the strength of the fuel region. 

With this condition the fuel element surface temperature will be 
less than the prescribed value of 850° F. The amount of damage that will 
result to the coolant from this temperature and the accompanying 
radiation is uncertain. However, for shipboard applications it does 
not appear to be a serious penalty. 

In the present day cruiser design for torpedo protection there is 


available a large amount of bunker fuel oil to partially fill the 
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layers of the torpedo protective system. Certain layers of the torpedo 
protective system are required to be liquid loaded. This system could 
be used for the storage of the organic moderator coolant. The unradiated 
coolant could be stored in the immer layers of the system and that 

which had been irradiated would be stored in the outer layers to protect 
the crew from any induced radioactivity. The long chain hydrocarbons 
would tend to settle to the lower parts of the tanks where they could 

be removed during periods of overhaul. With the installation of a 
purification Syeven du’ process of separation might be handled con- 
timiously aboard ship in the same manner that lubrication oil is 


purified and recycled. 
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Nomenclature 


A 
f 


2 
flow total area, (ft ) 
flow area per plate case n, (ft plate) 
heat transfer area, (£t7plate) 
average specific heat, (Btu/lb °F) 
hydraulic diameter, (inches) 
fraction of uranium in fuel region 
2 
mass flux, (1b/sec-ft ) 
: ' 2 
film coefficient, (Btu/hr-ft -°F) 
thermal conductivity of cladding, (Btu/hr-ft-°F) 
thermal conductivity of cladding, (Btu/hr-ft-°F ) 
fuel length, (inches) 
number of fuel plates 
total heat generation rate, (Btu/ft>) 


) 


heat generation rate per unit volume of fuel, (Btu/hr-ft 
fuel cladding thickness, (inches) 

coolant outlet temperature, (°F) 

fuel thickness, (inches) 


spacing between plates, (inches) 


bulk temperature of coolant, (°F) 
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Nomenclature (Contimed ) 


AT 


“ol 


fuel centerline temperature, (°F) 
fuel-clad interface temperature, (°F) 
fuel element surface temperature, (°F) 

i 2 
overall heat transfer coefficient, (Btu/hr-ft -°F) 
coolant speed, (ft/sec) 
mass rate of flow, (1b/hr) 
fuel width, (inches) 
overall temperature difference across reactor, (°F) 

| 3 
average density, (1b/ft" ) 


film temperature drop, (°F) 
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Basic Reactor Physics 


Criticality calculations were made using a one group reflected 
homogenous model with case 3, narrow fuel elements. Based on heat 
transfer considerations, inlet and outlet coolant temperature were 
set at 600° and 700° F, respectively, with a coolant speed of 15 ft/sec 

The assumptions in this calculation were 

le The structural members inside the reactor were of 
type 3:7 stainless steel and composed ten percent of 
the total reactor weight, 

2. A fast fission factor of unity, 

3. The uranium of fuel is enriched to 90% in yo? 

h. A right circulay cylindrical flux distribution of a cosine 
axially and a Bessel function radially, 

Se A beryllium reflector of 236 cm placed around the reactor. 

The reactor was treated as a homogeneous medium. By a one group 
calculation the volumes of each constituent of the reactor were cal~ 
culated from the dimensions given previously in the heat transfer sec- 
tion. These volumes were then converted to-weights and the volume 
fraction of the total reactor as shown in Table I. From the cross 
[5] 


sections given in BNL~325 effective cross sections were calculated 
for the operating temperature of the material. With these cfoss 
sections the reactor parameters as shown in Table II were calculated. 
These give an effective mitiplication factor of 0.89). This value is 
sufficiently close to unity that final adjustments are made by the 
computing machine. 

The primary purpose of the calculation is to obtain an approximate 


2h 





value of the critical radius which is needed as an input to the 


computer solution. 


Tahle I 


Reactor Materials 


Volume Mass 
(om? ) (Kgm) 
Fuel 6.95 x10" 498, 
U 6905 
2r 2565 
Clad ( Zr) 15.01 x10" 966 
Stainless Steel 1.99 x 10" 160 
Cr (18% by weight) 28.1 
Ni (11% by weight) 17.35 
Fe (71*by weight) 112, 
Coolant 26.3 x 10" 21965 


(73-5 Diphenyl oxide by yeight ae gil gor Pa «) and 


26.5% Diphenyl by weight (6 sledol 


C (06.91 by weight) 
H (06.03 by weight) 


O (86.59 by weight) 
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Volume fraction 


0.0068 
0.108 


0.007. 
0.00) 
0.027 


0.4703 
0.03292 


06377h 


Table II 


Reactor Parameters 





5 -# Beryllium shield thickness « 2306 CMe 
h Las 1 
om oo Vv — —S1 = * 
a. z, 2 (1 Ho), fh 5 Dy 0.2833) cm 
' x rr m 
Di 3 =- = 10176h4 em 
te we 
= i, +B = t-- 4 
d 0.71 P| te ip” 3m 39.599 cm 
| : ai 
5 z »F V = 0.1208 cm 
ay f 
i 
1° is . 97387 tn’ 
a 
2h \ 415 
p(E) = exp Nog (369) - am x 10-4 = 0.726 
2. NA 6 
an 
1? A) = Si 1.9692 
a{ total ) 
k = HP S 1244296 
Ea. § 
1 “ijt 
_ = 0017598 
“ls 
D : 2 
= inj — = 343 cm 


2 2208 ~3 
B (ee 7 = ete 1g332 x 10° cm 


B = albedo of Beryllium 
-Ber 


k@ 
Kare 14+B2L2 0.89h 
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Application of the Flux Equations 


The S approximation is a method of differencing ih which the 
independent variable is quantized into n intervals, and intermediate 
values of a function of the variable are approximated by a straight 
line between points at the boundaries of the interwale This method 
of numerical integration was developed by Bengt G. Ssvtgdal's 

The vector flux at a point within a reactor can be obtained by 
the use of the 5, method. The equations necessary are derived in 
Appendix A for two cases of neutron velocitys 

Case A - The component of the vector flux in the axial direction 
is greater than its component in the radial direction. 

Case B - The component of the vector flux in the radial direction 
is greater than that in the axial direction. 

These two cases are illustrated in Figure I using a grid which is 
used to difference the transport equation in the r and z coordinates, 


r being the radial coordinate and z the axial coordinate. 


Figure I 


The two cases of flux ata grid point 





Z¢ 
Direction of 
neutron motion 





Cass B 


Direction of 
neutron motion 
Case A 





The general flux equation for case A iss: 


[». | anus? (‘et ,AgLe gobt Ab 





LULU 4 7 z coal jm 2 
J ym igh. Cc. 
>,| td (A >) ? 


where h is a value which is one greater or less than j, depending 
upon the direction of the radial component of vector flux, and n is 
a value which is one greater or one less than m, depending upon the 
direction of the axial component of vector flux. Definitions of the 
other quantities can be found in the nomenclature part of this 


sectione 


For Case B the flux at a point is defined by the following 


equations 


p talenah (gs Naat Mes AA) 
Gl ( ¢.* 2 " a 


‘i - 








ni 


(2) 


where A for = (1) and (2) is a bys: ) 
\s Cc 
. a “4 is An “3 \Mi-1,f 
A=\|p,)-5 (, a, "§, “jot bls 2 (3 lhe 
A ie : A 
. p,| lede 3 - | ami,t , °. E-| A. ay i 


hh 2 jgm h,n 


Hack lene a. A e? G3) 
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It should be pointed out that the total macroscopic cross 
section, 2, is evaluated using the material properties at the 
center of the position grid of Figure I, but are designated by the 
grid values at the corner diagonally opposite the point j,m, 
iee.3 the point h,n. Further, rs is actually the mean radius and is 


equal to j oh . 


2 


The source at a point is defined by the following equations 
k k kop, 
4 -“ Tee ae $ Ak a 
jm kel] f “jym kel 4 ji 
where the scalar flux e J is equal to 
3 
A gist 
H 3 

f=0 fiso i oe 
The fission cross section and transfer probabilities, C, are 
evaluated using the material properties at the grid point at which 
the flux is being calculated. The data necessary for the solution 


of the flux at a point is summarized in Tables I. and Il 


a the S 9 approximation 

: c d e P 

as = = _— Las =i i 

0 “1 -1 0 = 1 40.3183 
1 0 -1/3 14/3 -2/3 2 #03613 
2 1 2/3 ),/3 1/3 2 +0. 3183 
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Data for Gauss quadrative approximation 


Table ITI.. 


ft - = os “t 
0 ~0.8611363 1.691625 #1.967705h +063),785))8 
i -0.3399801 0.319728) #1063310 +0.6521))52 
2 +0.3399801 +0. 319728), +1.0633)0 +0.6521)52 
3 +0,.8611363 #1.69))625 #1.967705) +0.31:78518 


Not all pairs of values for i and f are r@levant for each of the 
cases A and B as a consequence of the defining inequalities for these 
cases. The pairs of relevant values are discussed in comection with 
Figure Il, A general point may be defined to be any grid point of the 
reactor other than a point at the center, on the centerline or in the 
midplane of the reactor. A point at these other places is called an 
origin point, a centerline point or a midplane point, respectively. 
For Case A for which|b | <|h,| and in the case of the first quadrant 
flux, there is a further requirement (Table I, Appendix A) that b, 
and H have negative signs. In this quadrant there are three Case A 
flux equations, one for each of the pairs (i=0, f=0), (i=1,f=0), and 
(i=1,f=1). Similarly, for Case B, the requirenent|b | >|h, |leadg to 
only one relevant combination. 

The applicable combinations and those for the remaining three 


quadrants are summarized in Table ITZ. 
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Figure Il 
Flux Configuration 


for a particular energy group 
at a general grid point 


v- 


m+} 





Table III 


Allowable combinations of i and f 


Quadrant Case A Case B 


i=2, f ‘=2 





The configuration of the flux at a grid point which results from the 
allowable combinations of i and f are shown in Figure II. If the flux 
is assumed to be distributed inXe1 energy groups there will be 12 (2 #21) 
flux paniliiwi adtoeieted with every grid point. For the case considered, 
there are three energy groups and, therefore, thirty-six flux equations 
for each point. 

The following boundary conditions are assumeas!@? 

1. Continuity of flux exists at the centerline and at the 
midplane of the reactor, and further the flux is symmetrical 
about these boundaries. This permits the calculation of 
the flux in only one quarter of the reactor. Figure III 


shows the flux distribution at these various special points. 
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Figure III 


S, approximation to the vector flux distribution 
at the boundary points in a reactor 


2 





2. A second boundary condition is that at the outside 
boundaries of the reactor there exists no inward flux. 
In order to satisfy this requirements, an extra set of grid points is 
placed just outside these boundaries, and at these points the vector 
flux whose argument describes inwardly moving neutrons is assumed to be 
zerde Note, however, that leakage does exist at these fictitious 


points and at the outside boundary points of the reactor. 
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Nomenclature 


b, 
i 
¢. Values obtained by differencing: and: integration of 
a. the transport equation over 4) pg and defined by 
ms 

e. Equation (7a to 7d) of Appendix A 

ee .Yaluas depending upon the quantized values of p and 

h i defined by Equation(7e and 7f)of Appendix A, 

it, Values obtained by the use of the Gauss quadrature 
approximation to sum the flux over p spangmee . 
defined by Equation (13) Appendix A. 

us Values obtained by the use of the 5, approximation 
to sum the flux over ')/ space, defined by Equation (12) 
Appendix A. 

r mean radius between two adjacent grid points 

S Neutron source 

4, Fraction of fission neutrons which are produced in a 
particular energy group ibe 

ov Number of neutrons produced per fission. 


3h 








Le 
rt 


A(r,z, ty Vv) 


The probability per wit length of neutron travel 
that a neutron is scattered in group k and that a neutron 


emerges in group 4 from the scattering event. 
Macroscopic fission cross section 
Total Macroscopic cross section 


neutron fiax which is classified according to speed, 
V, colatitude and azimuthal direction, cosines p and ¥, 
and position r and Z. 


spacing of grid in the radial direction. 
spacing of grid in the axial direction. 


cosine of the angle © between the axis and the 
direction in which neutrons move. See Figure I 


Appendix A, 
Cosine of the angle We See Figure I Appendix A. 
See defining equation, page A-12 Appendix A, 


See defining equation, page A-17 Appendix A. 
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Nomenclature 
Superscripts and subscripts 


i index labeling the quantized values of ]] ° 
x index labeling the quantized values of pt. 
3 grid position index in the radial direction. 
, \ 
m grid position index in the axial direction. 
L neutron speed index. 
k ~~ neutron speed index of group in which neutron is 
scattered. 
n axial prid position diagonally opposite grid point 


at which flux is being calculated. 


h radial grid positions 
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Ack 
Computer Introduction 12] 


A computer program is a sequence of logical instructions telling 
the computer what mathematical operations are to be performed on the 
data. Therefore, there are two groups of inputs, the data and the 
instructions. The data fed into the computer are stored in the fast, 
core memory if there is sufficient roome When the space is not 
adequate, some of the data are stored in magnetic drums. the only 

_ difference in the two storage processes is that the data in the 
drum takes longer to locate and bring into the core memory. Each 
number or piece of data is stored in a register or cell the location of 
which is designated by a number called the address of that register. 
The instructions are stored in sequence and are performed in the order 
in which they are stored mless the program directs otherwise. No 
unused registers are permitted to interrupt the storage as these will 
stop the computer. 

The instructions given to the computer consist of two parts, the 
operation and the address sections. The operation section is first 
and uses letters to tell the computer the mathematical operation to 
be performed. The address section gives the address of the register 
containing the data on which the operation is to be performed. The 
instructions recognized by Whirlwind I are listed in the glossary. 

_ For Whiflwind I there are many useful devices that make the 
programming much simpler. A few of these devices used in this 
program are floating addresses, preset parameters, and subroutinese 

Probably the most useful device is the floating address. This 

floating address or "flad" is used to eliminate the necessity of 


referring to the absolute register address in the core memory when 
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writing the program. The absolute register address of any part of 
the program must be known primarily so that the initial or final. 
instruction of data is not placed in a register containing part of 
the subroutines fixed in the memory of the computer. 

In the program key instructions on which mathematical operations 
are to be performed or which will be referred to later in the program 
are given floating addressese These addresses consist of a letter and 
a number, for example al, a2, etce Any letter except o and ell 
may be used and the sum of all the numbers following these letters may 
not exceed 255. When the program is read into the computer, it will 
assign absolute addresses or register numbers to these flads. A 
table of the floating addresses is produced by the computer with the 
solution unless it is specifically suppressed in the computer. This 
table gives the absolute address of every assigned floating address 
and is particularly helpful in trouble shooting a faulty programe 

Another useful device is the preset parameter. These are 
numbers other than data that are used in the program. Preset 
parameters are frequently used to designate the number of times that 
a calculation is to be performed. Preset parameters mst be assigned 
a value before they are used, and they will retain this value until it 
is specifically changed. Preset parameters are designated by two letters 
and a number. The first letter must be a py u or 2 and the second 
letter and the number may be anything, except the letters o and el may 
not be usedo ; 

A third useful device is the subroutinge Its usefulness comes 
from the fact that it is a complete entity in itself and may be construce- 
ted and trouble shot independently. For the build up of a complete 
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program it is possible to start with one of the most basie mathe- 
matical operational sequences that will be performed repeatedly, 
program it, perform the hand calculations of the sequence and 
then place the routine on the computer. Modification of the 
subroutine or correction of the hand calculations mst continue 
until the results of both are in agreement. 

In arranging these subroutines in the main program extreme care 
must be used. The main program must be interrupted to call the sub- 
routine into use and after the subroutine has performed its function 
it mst direct the computer back to the proper place in the main programe 


An example of a subroutine is as follows 


Main program 
isu p 10 
tl, ispx take the next instruction from register 
x ) and contime from there 
ime a 5 forms product of (-m2) and (a5) 
i STOP end of program 
xh, itax § transfer into the address section of the 
instruction in register x5, the address 
that is one more than the address tl. 
ica m 2 ; 
places m2 ina 7 
its a 7 
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ics m 2 
places =m 2 ina 8 


its a 8 
x5, isp 0 returns to main program one register 
after t 1g ie., register containing 
imr a 5. 


Since the numbers used in this program were over four decimal digits 
in length, floating point arithmetic was required. In this system two 
registers are required to express a number. Each of these registers is 
capable of containing 16 bits of binary number information. The first 
bit of the first register expresses the sign of the number and the 
first bit of the second register expresses the sign of the exponent of 
20 Six bits are required to locate the binary point of the mumber 
leaving 2) bits to express the actual number. In certain stages of the 
program it is desirable to handle single length register mmbers. The 
instructions for the handling of these numbers are similar to those for 
double length registers except that the i does not precede the instruc- 
tion. When operating on double length registers the computer is said 
to be in the interpretive mode, when on single length registers it is in 
the Whirlwind mode. To enter the interpretive mode from the Whirlwind 
mode it is necessary to use the instruction "IN" and to reverse the 


direction it is necessary to use the instruction "OUT". 
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GLOSSARY OF INSTRUCTIONS 
TERM DEFINITION 


ica x clear MRA and add contents of 
register X. 


ics x clear MRA and subtract contents 
of register x. | 


jad x add contents of register x to 
what is in MRA. 


isu x subtract contents of register x 
from what is in MRA, 


imr x multiply what is in MRA by the 
contents of register x and round 
off product to fifteen digits. 


idv x divide contents of MRA by the 
contents of register x. 


its x transfer the contents of the MRA 
to register x losing what was 
previously in register x. The 
MRA remains unchanged. 


iex x exchange the contents of the MRA 
with the contents of register Xe 


isp x take the next instruction from _ 
register x. Does not affect the 
contents of the MRA» An uncon- 
ditional instruction used to break 
the sequence of operationse 


~ (= « 
_ £* . 
> ed 
= 
= 
i 
=> 
>_> > 
=_ = 
7 
> 
> = 





Glossary of Instructions (Continued) 


ep x 


ise x 


icr x 


ict x 


ici x 


icd x 


DEFINITION 


if the contents of the MRA are 
negative, take the next instruc- 
tion from register x and continue 
from there. If the contents of 
the MRA are positive, take the 
next instruction in sequence. A 
conditional instruction used to 
break the sequence of operations. 


select counter mumber x. Without 
this instruction counter zero 
will be used whenever +C appears. 


cycle reset. Sets index register 
of counter to zero and criterion 
register to Xe 


cycle transfer. Increases con- 

tents of index register by onée 

If contents of index register greater 
or equal to contents of criterion 
register, set index register to 

zero and do next instruction 

in sequence. If contents of 

index register less than contents 

of criterion register, take next 
instruction from register Xe 


increase the contents of index 
register by number Xe 


decrease the contents of index 
register by number x. 





Glossary of Instructions (Contimed) 


TERM DEFINITION 
iti x transfer the right eleven digits 


of the index register into the 
right eleven digits of register Xe 


iat x add contents of index register to 
the contents of register x and 
store the result in the index 
register and register Xe 


ita x replace the address section of the 
instruction in register x with the 
address that is one more than the 
address of the register containing 
the last isp or icp if the contents 
of the MRA are negatives 


iLTOA record the contents of MRA on 
direct printer (typewriter). 


i MOA record the contents of MRA on 
magnetic tape for delayed printings 


iSOA record the contents of MRA on 
oscilloscope for photographingse 


iFOR this instruction provides an 
automatic device for obtaining a 
suitable layout of output data 
in colums, lines or blocks. 


Notes The xbove output instructions are usually followed by a series 
of letters and numbers that indicate the desired form and 
arrangement of the outpute 
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C ter Logic 


The program starts at 1 in the main program by entering the 
data, program and preset parameters. The initial calculations are 
with the reactor radius equal to the value placed in register rhe 
With this radius the program then shifts to the flux iteration and 
convergence test, 1. In this routine parameters are set to control 
the number of iterations allowed at one radius and to traverse the 
computer thru all points in the reactor. With a value of j and nm, 
defining a point in the reactor, the program shifts to Flux, 1.5, 
where the addresses on the drum of the required data are calculated, 
the data entered in the fast memory am the routine appropriate to 
the point is selected and control is shifted to that routine. 

The first calculation at any point is the scalar flux for all 
three groups and the second is the input from fission and scattering 
from an upper energy group. The first point to be calculated in each 
eyele is the origin. This is the point of normalization of the 
inputs to the fast group. The normalization constant, 0 , is defined 
as the reeiprocal of the fast group input at the origin. With this 
constant all other inputs (all groups and all points) are normalized. 

As explained in the previous section of the procedure certain 
quadrants of the flux are required for each point depending on its 
locations Prior te sequencing to another point all fluxes appropriate 
are calculated, 





The fast source input is then compared to the previous fast 
source input. If the new value is not within a set specified fractional 
value of the previous value a negative quantity is introduced into s 10. 

When the entire grid has traversed a check is made to see if the 
maximm allowable attempts at source convergence have been made. If this 
maximum muber has been made the computer will stop. If not a check 
will be made of switch s 10, As indicated in rogtine h, another eyele 
will be made thru the reactor if # 10 is negative, and if s 10 is positive 
control will return to the main program to repeat the above procedure 
based on the radius placed in yr 3. 

When a gamma is calculated for this radius, control is again 
returned to the main program but at a different point. All subsequent 
returns are at this point. A check is made to determine if the 
maximgm allowable molbér of iterations of radii have been attempted. 

A test for criticality is made to see if gamna is within a 
specified tolerance of mity. If it is not, a new radius is found by 
extrapolation or interpolation based on the most recent values of R and 

¥ and the best previous value of R and J. In this determination of 
radius, the gamma that is closer to unity and the associated radius are 
preserved for use in later extrapolation. The previously outlined cal- 
culations are repeated with this radius. If X lies within the specified 
tolerance of wmity control remains in the main redting and the print-cut 


routine, which has been stored in the drum, is entered into the fast 


menorye 
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The computer is then sequenced thru the grid calculating and 
printing out scalar fluxes at each point for each groupe 

If the flux pattern is not as desired the properties of the regions 
are changed and a new tape of data cut. The program is then run again to 


determine the flux pattern. 


L6 








{Co 








Enter preset 
parameters 


deta, Roan R 
bh heii a 


ip 


Flux iteration and 
convergence test 


R=R 
k Flux iteration and 
convergence test 


5  R iteration 
icheck, n'3 N® 


Test for 
criticality 






out routine 
in MRA 





Li7 


Main routine 







Flux iteration and h 
convergence test 

cea ee Teil if 

Enter from main routine a) 

or radius extrapolation (9) 


4.3 |Set source convergence 


switch (s10) 










Grid 
incompletely 
traversed 





conver gence 
switeh ~ 





Exit to 

R=R_ (5) on first cycle, 
to iteration ad 
convergence test (1,) 

on subsequent cycles 
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aux hes 
Enter from Set J,M he 
grid (ek) or Grid 
Traverse (1.6) 







alculate addresses o 
required data in drum 
d enter in MRA 









a Throw 
~~ | source convergence 
switch (set 510), 
if required 





Exit to 
Grid traverse (1.6) 
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Grid Point Sensing 

Enter from Address calculation and 

data entry (1501) or calculate drum 
address of required data and enter 

data (8.1) 


Centerline 
Program Program 
(4503) (505) 
or or 
tae Centerline 
Scalar 
Scalar 
Flux Finx 
(453.1) (48502) 
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Midplane 
Program 
(450) 


or 
Midplane 
Scalar 
Flux 
(45401) 


Interior 
Program 
(4506) 
or 
interior 
Scalar 
Flux 
(45601) 





Origin program 4503 
Enter from Grid Point 

Sensing (5.2), origin 

exit 









45301 
Print (8.2) BS a 
\y-2| s Ey 
45302 
45303 
530k 
45305 
Jom 
ist Quadrant 
oh ee 


Throw source cone 
vergence switch 


(4507) 





Midplane program Sols 
Enter from Grid Point 


Sensing (5.2), midplane 
exit 











USliel 
——>- _ Print (8.2) if 
|¥- s ey 
115302 
LS3ok 
45305 
Sho 
Exit to 


Throw source cone 
vergence switch 
(4507) 


52 


‘4 





Centerline program L505 
Enter from Grid Point 


Sensing (i502), centerline 
exit 


Print (8.2) if 





Yul s ey 
hth Quadrant 
| Exit to 
Throw source cone 


vergence switsh 
(507) 


oD 





Interior point program 
Enter from Grid Point 


sensing (15.2), general 
point exit 


4506 


45601 [9 
jgm —> Print (8.2) if 
ali £, 5#0,mf0 \y=2| ee 








A gtot 
Jom 
ist Quadrant 


oho | gif 
om 
2nd Quadrant 






4562} 2g ra 
aa 
3rd Quadrant 


gist 


155.2 


Jom 
hth Quadrant 


Exit to 

hrow sources 
convergence 
switch (115.7) 





5h 





Tiwow source convergence switch a 
Enter from origin, midplane, ° 


centerline or interior point 
programs (15.3 to 5.6) 





I 


Exit to 


Place flux in drum 


(4508) 


op 





Grid Traverse ate 
Enter from Flux (1.5) 





Exit to 


to Exit to 

x iteration check (1.7) Fiux (5) or 

or Stop Drum address 
calculation 
(851) 





Print out routine 8 
Enter from place print 


out routine in MRA (7) 


hel Set Jy M 
Grid 


Bel Calculate address 


of required data and 
enter data ‘ 


502 Grid point 
Sensing 


ye 
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RESULTS 


Numerous runs of subroutines and of the entire program were made 
on Whirlwind with the program contained in Appendix B. This program 
was adapted for a reactor radius of 75 and 95 centimeters and a height 
of 168 centimeters. The grid placed on the reactor was a 6 by 6 grid 
giving for the first trial radius a grid size°15 cm by 3) cme 

The final run had programmed a first trial radius of 75 emg “— 
maximum allowable number of iterations with a constant radius was 
twenty and the tolerance of convergence of the fast source input was 
three percent. The computer stopped, as programmed, upon completion 
of twenty iterations at a set radius when the fast neutron input at. 
all points was not with the tolerance specified. The computer running 
time per iteration is two mimtes. The computed values of scalar 
flux at each point and energy deeoup were printed out for each 
iteration. The normalizing constant decreased from 0627 for the 
first iteration to 0.01 for the twentieth. Some values of the scalar 
fiux at the following point were negative 

a) interface of core and reflector, 
b) reflector, 


~ 


c) outside boundary of reflector. 
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DISCUSSION OF RESULTS 


The negative values computed for scalar fluxes are known to be 
unrealistic. The possible causes of this error are 
1) improper programming of the theoretically derived 
flux equation, 


| 2) grid size being excessive. 
tthe theoretical background for the flux equations is sound, and 
these equations have been programmed for other geometrical deen ied te 
on other computers and give correct ‘flux distributions = The 
need for a complete hand calculation is well realized, but due to lack 
of time these hand calculations were not completed fully. Therefore, 
there is no reason to assume or expect that these equations have been 
programed correctly. To accomplish the required hand calculations, 
computer verification of the hand calculations, and program trouble 
shooting it is estimated that another three to four weeks full time 
would be requirede 

For accuracy in evaluation\of neutron transport equations, the 
grid spacing should be Less than or approximately equal to the mean 
free path of the nevtrone If this condition does not prevail, the 
accuracy and reliability of the calculations decrease rapidly. 

A two fold increase in the number of grid points will increase the 
time of computation by a factor of eight. sey increase in the time 
required for one iteration thru the reactor is not allowable since ndw 


two mimutes are required. 
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The code developed is in the 'interpretative' or ‘floating 
point’ mode for Whirlwind, for the reason that this mode is simple 
to learn and to trouble shoot. Whirlwind has also available a 
Whirlwind’ or 'fixed point! mode in which the speed is increased 
by a factor of roughly 35. Use of the fixed point mode would 
make the two dimensional S, calculation on the Whirlwind computer 


practical. 





CONCLUSIONS 


1) The code as presently programed is incorrect. 

2) The grid spacing as used is much too coarse. This 
difficulty is not inherent in the coded programe 

3) The mode of computer operation is unsuitable for an 
iterative process such as this because it results in 


a gode that progresses too slowly. 
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RECOMMENDATIONS 


1) Prior to the use of any computer a formal or semteforma] 
course in programming should be completed. 

2) A grid spacing comparable to the mean free path of the 
neutron should be selected. 

3) The fastest possible mode of calculation on the computer 
available should be usede 

h) For a coded program of this difficulty a faster machine 


with a much larger fast memory is very desirable. 
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Appendix A 


Development of Equations Using 55 Approximation to the 


Boltzman Equation for Cylindrical Geometry 
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Equations will be developed for the time independent case using 
cylindrical geometry. 

The vector flux at a point in space is a function of two variables 3 
one describes distance along the vertical axis of the cylinder, and the 
other describes radial distances. The axial variable will be designated 


as wp and the radial variable as ne From Figure I it can be seen thats 


n = cos Y ; 
Figure I 
u = cos © and , - 
(a sa omponents of neutron flux 
sind= V1 = at a space point. 


sin = Va - ne 

aN 
and where £2 represents the unit 
vector in which direction neutrons 


travel at a point. Therefore, 


ag- 328, (2 2855 ; 22 


rar Or dv ere 





a ~ iu 
and 2 = ji sin 9 eos 
r 
“an 


3 5% hone 
@ sin 9 inp =f i Cos Oo 


Therefora 
: | e 8 e Pe a : 
A .g = sino cos? 2. in osind 29, 06628 | (1) 
ar r oP oz 


If we substitute the expression for al and « which were derived from 


Figure I in equation 1, we find 
Vea (1 - @-4°) | o- ¢ 
E's Vg = ni bt 22, l= ze 4 2h (2) 


The quantity pe % represents the leakage of neutrons out of an 


element of volume per unit time in the direction A per unit solid angle. 


Keo] 





In the time independent case this leakage plus the rate of loss of neutrons, 
>» @, per unit time must equal the source strength, S. In problems 
in which cross sections vary with energy and neutrons alter their speeds, 


[6] 


the multigroup approximation of the Boltzmann equation is introduced: 


wk C yess r 5 a’ 
28 , ken 28, 2b ab - ob, (3) 
Or r an View a2 Yiu" Tap *” 


} 


where As is the total cross section within the energy group apt and 
ads is the neutron source term which results from fission and energy 
degradation of the total vector flux at the point. Slowing down is 
incorporated into the source terme 

Equation 3 is identical in form to the time dependent multigroup 
equation for spherical geometry in which the z term is replaced by a 
term representing the rate of change of flux with respect to time. 

The = approximation is a method of evaluating integrals over a 
variable in which the variable is quantized into n intervals and 
intermediate values of a function of the variable are approximated by 
a straight line between points at the boundaries of the intervals. 


We shail consider ¥( Space as the first variable which will be inte- 





grated. The flux at any point Figure IT 
(P9Z59 ) in Figure II can be Th Mesh 
Bir,z,m%) Actual | 






represented by the following flux 





equation where i represents 


the interval number. 


Aa2 





l li fe] ie A 
= ae O(r,25 ) SS B(r525 ; ) 
Ge aa Li a Lied (4) 


inl 
“For the S, approximation i = +1, 0, and -l. If we denote Nh ad Nant 


P(r 52> l ) = 


as A and par ) as g we obtain 
( eas) 2 
i - (r a g ” xe ~ AL g 74. 


a i 


g = 


If we substitute this in equation 3, we find 


is a shacn ue a ag boa e tua! ist 


Ay or 24 On 
al : ore -n) me re Wn - (ew an — a4 x: 
On ” Vien 4. 2 a, dz 
2 Lash ZL io] 
= | Ay = sy) g zn on) a (3) 
{iepe i. 2 \ low | - a. \ Tous 


Since the source, §, is isotropic, integrating this expression over if) 


from n fel een 2 we obtain n equations in the form 
2 1 
LA h 43 io] 
Pw o, 225 « Q ew, EE. ~~ 
ir 0 es 
2 Aol , 
on2B, Ht. es. (6) 


where 


; 3 (7a) 


ce = (3 ; m 
__” AG ae eee (7b) 


Ax} 





(7c) (7e) 


o 
a oe ee 


¢ = 2 (7d) g.= : (7f) 


A different method of finite differencing is used for the other 
directional variable, A 9 Since use of an expression analogous to ()) 
for in the integration of the equation (5) with respect to p would 
result in an extremely complex expression. Instead, the equation (5) 


is directly quantized with respect to p giving the results| 


ey ee ce «yee 
ee 


iar “i r Me OC 2 
(8) 
argirlof d i-l,f Ralf 9 pg. 
a Se on Balen the sa ae 
i@ge . f @2 is = 
We may now proceed to evaluate the source integral at a point. 
For isotropic scattering arm three energy groups of neutrons 
p k., Ak 
s(rj2) =p 2 2, Broz) + BT Mroz)s 
(9) 


where 
+) +1 
| aN 


"soas) = 5 fa { Sevtomin pm 


As in equation (),) B(r 5%) is approximated by 





u) + 
4. 4, 


Me. eyes) 6 EL 


P(r Zo tty )e g(r ,2 5h), 


A=}, 





and 2 Ie 


i (c! 
k 1 ] d d 
P(rozo) = L. | / aria (core | 


Ot ij=1 





Upon integration 





; pi pa oe 
= ‘Le o \} Leo + sin 
P(r 5% 51) ; > |, | Nana il 1G 1 eel 
4.20 
=] 
_ ol Sloe g (11) 
=! gaeee a "9 vm A Ns “? a 
i 
in da ; \ 1 
a sir o- 
i j it iol Li (Ben : 
If values of ns = +], O and -1 (S, approximation) are substituted in 
the above expression, we find that 
k a 4A 
Blro%yu) & On e mA B (r 5251) 9 (12) 


where P. = 0.3183, an = 0,363, and P, = 0.3183 . 


By use of the Gauss quadratdive method 8! the remaining integral 


can be evaluated in the form of a sum of the products of coefficient 


times the vector flux where H values are the nt2 positive roots of 


Pp = Oo 
ae @ 


A-s 





The table below lists the values of and the coefficient a used 


f U H 

f f 
O ~0.8611363 053k:7858 
al =0. 3399810 0,6521))52 
2 +0. 3399810 0.652152 
3 +0,8611363 0.378518 


Therefore, since H = e and EF = H, the scalar flux at a point is 


Bez) 3 H $ P ‘fF (42), 


f=0 f i=0 i (13) 
and the source term equals 
4 k k Ak k 
S= 7 3 > + 3 x ZB )o - ; 
p ke] e g (r,z) k=l / Kr,z Die (14) 


where 
" is the propability that a fission neutron will have an energy % and 
Ak the probability that a neutron having an energy k will be 
scattered into the energy group Z , k> fo 
The scattering probability requires further explanation. We define 
Lb Lk as the probability that a neutron which is scattered in the ie 
group will land in the Gn group, k as the energy group in which the 
scattering occurs (1 for fast, 2 for epithermal and 3 for thermal) and 


ZL as the energy group into which the neutron is scattered. The 


energy range covered by scattered neutrons is E-aE, and the proba- 





AE 
bility of landing in the interval AE is 


S E 
E-ak Let E be the energy 


of the neutron before scatterings 


B <E4E. 
ae 


A= 





gg 


The neutron flux per unit energy interval is JE ° 


d 
Then id - dE = (E) is the mumber of neutrons scattered per unit volume 


d& Ss [9] 
and per unit time; but 


gig 
dE SE F 
s 
d q dk 
Therefore Y a 4 (BE) “oS = af 
9 4, GE & (8) ; (15) 
E -E 
Vy £4 or ; : f 
~ E(ea) is the probability a neutron will land in the th groupe 
Lk 


‘ris the number of neutrons scattered in the k th group that land 
in the £ th groupe Since the absorption, sources and leakage in a 
large thermal reactor are small, the slowing down density q, is rather 
constant over an energy group. This fact is used to compute the num- 
ber of neutrons scattered from the group k to group Pe f10] 


The equations for each case are listed below. 


Cage Equation 








KS 
lo E, > aE aT —— a Ya - Se) (16a) 
Jel ok (Lea) ine) . Ba Es =< E 
KS E =i E 
26 Ey>aek Ak Ss | ae ee 
Te) ee TR 
3 1 Sco? k¢1 k 
k L+l k+l 
ak 
= @ Te < (16b ) 
£+1 


A=7 
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We may now proceed to evaluate equation ] by integrating over the 
r and z variables. However, before doing so we must consider the direction 
in which integration should be carried out. Itis always best to integrate 
in the direction of neutron flow since any errors produced in the process 


of integration will be diminished as we proceed further in this direction 
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If integration is carried out in a direction opposite to that of neutron 
flow, any errors produced will be magnified. 7] 
We mist now difference our equations with respect to r and Zo 


We can arrange terms in equation (8) in order to facilitate the integra- 


tion. Thus equation (8) becomes 
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We mist now further divide the results of the integration into two 
Separate cases?s 
Case A = The component of velocity in the r direction is less than 
the component of velocity in the z direction. 
Case B ~ The component of velocity in the r direction is greater 
than the component of velocity in the 2 directiono 
If we consider Case A, and difference r and z ag illustrated in the figure 
below, we obtain four distinct subcases. These will be classified in 


accordance with the direction of 


Figure III 
Position Mesh travel being determined by the signs 
nl of Ly and He A positive sign 
indicates a direction of travel away 
from the center of the core. The 
va various cases are listed in Table I. 
me 
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Table I 
In all A cases fb, | <(h, | 
In all B cases|b | >{h, | 
i ft 


Subcase a fe 
Al, Bl - o 
A2, B2 + ~ 
A3, B3 + + 
Aly Bh - + 


Considering subcase A-3 and integrating equation 5 and knowing 
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Figure IV 


| and 
Position Mesh for Case A 
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we obtain the following equation 
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Where ie a is the total cross section of the material at the center of 
j= 


the grid in figure lg i.eo, the material to the left and below the j,m 


= er. 
point and, | is rj tel 
a ° 


The radial direction of integration is determined by the sign of rl F 
1o@o, by the sign of the ae or .° Consequently, the signs oe 9 by» 
d.» and 4, all change together. Likewise, the direction of integration 
axially is determined by the sign of Hs Leee, by the sign of a? 
Consequently the signs of un, h., and " all change together. If we 
use absolute magnitude signs for these quantities all Case A equations 
will be identical insofar as the signs of various terms are concernedo 


Similarly, all Case B equations will be identical. 
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If we mow add the following expression to both sides of equation (18) 
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And finally transposing terms in equation (18b) we obtain an expression 


for the vector flux at the point j,m. 





By use of similar procedures for developing the equations for the 
three remaining A subcases we obtain the three equations on the 


following pagese 
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For subcase A-}; 
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For subcase A-1l 
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A similar procedure is used in deriving the equations for the 


B cases$ except that as shown 


Figure V 
Position Mesh 
Case B 


in Figure 5 the r and z com- 


ponents of the flux are different 














because b, phe The method 


of deriving the flux equation 


by 
for case B=3 is outlined on 
the following page. The 
equations for the other B cases 
are listed also. 
m=1 
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To both sides of equation (22a) we add 
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where . is defined by equation (19b) 


For case B=] 
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where A, is defined by equation (20b) 
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For case Be2 
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where A, is defined by equation (21b). 

The vector fluxes defined by equations (18d), (19a), (20a), (21a), 
(22c), (23), (2h), (25), are valid for all values of u with the excep- 
tions of the end points of our w space where 1 = +1. At these points 
trouble is experienced with singular evoke newereey a modified equation 
may be derived by letting the terms scimatios © 4 all equations 
following equation (1) be replaced by zero, since for » = + 1 the 
term ia in equation (1) vanishes. This special equation was not 
needed since values of were chosen which avoided these singular 
points. 

The procedures described above result in n equations. One more 
equation is needed and is provided by setting n = <1 in equation (ye 
The development, above may be used when the values of b, c, 9., a4, 
and @ used above are replaced by -1, 0, Z( Nw) O, and 1, respectively. 
These values reduce equation (6) to that found for the special case. 


Thus equation (6c) becomes with i = 0 
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Nomenclature 


et 


Values obtained by differencing and integration of the transport 
equation over ” >» and defined by Equations 7a to 7d of Appendix 
A. 


Values depending upon the quantized values of 4 , and defined 


by Equations 7e and 7f of Appendix A. 


a unit vector in the direction indicated by a subscripte 


Fraction of fission neutrons which are produced in a particular 


energy group Xe 
The radial distance to a point. 


The mean radius between two adjacent grid points. 


The axial distance. 
Neutron energy frem a point. See Figure il, Appendix A. 


Values obtained by the use of the Gauss quadrature approximation 


to sum the flux over space, defined by Equation 13, Appendix Ae 


Values obtained by the use of the S . approximation to sum the 


flux over space, defined by Equation 12, Appendix A. 


Ratio of the minimum neutron energy after scattering to its 


energy before scatterings 
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A Grid spacing of a mesh, 
4] Cosine of the angle WV 


e Angle between the axis and the direction in which neutrons 


travel. See Figure I, Appendix A, 
UL Cosine of the angle 9. 


th 
HO, The probability that a neutron which is scattered in the k 


th 
group will land in the @ = group. 


a) Number of neutrons produced per fission. 
D> Total macroscopic cross section. 
Fission macroscopic cross section. 

Qk 


V The probability per unit length of neutron fravel that a 


neutron is scattered in group k. 


% (r,z, Neutron flux which is classified according to speed, v, 


a (i Y) olatitude and azimuthal direction cosines u = and 


position r and 2. 


Y Colatitude angle denoting neutron travel with respect to the 


radius. See Figure I, Appendix A. 


OO See defining equation Page A-12, Appendix A. 


a 


See defining equation Page A-17, Appendix A. 


@é 


b> 


Direction of wit vector in which neutrons travel, 
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Superscripts and Subscripts 


i Index labeling quantized values of ie 

f Index labeling quantized values of i’ 

3 grid position index in the radial ities 

m grid position index in the axial ancl 

Q neutron speed index. 

k neutron speed index Be group in which neutron is ee 
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Prior to the use of this program the following preset parameters 
must be defined and/or calculateds 
pzl = number of grid points in a radial direction less one 
pz2 = mmber of grid points in an axial direction less one 
p23 = number of iterations, N, allowed at a set radius, in 
attempting to achieve a convergence of source patterne.-. ; 
p2l,; = number of adjustments,iN', of radius.allowed in 
attempting to achieve erties itty. 
pre = drum address of cross section 
pfl + 92(pz1 + 2) (pzl + 2) 
pf3 = 92(pz2 + 2) 
pfl = 6(pz2 + 1) 
The preset parameters listed below were used in the development of 
the program and have the values as indicated 
pf1l = 5000s drum address of first flux at origin pointe 
pf5 = 703 address in fast memory at which the print-out 


routine will be placed when needed 


\ 


pf10 = 26503 address of the print-out routine in the drum 


pf15 = 3003 length of the print-out routine 


pf13 = 11,703 address in the fast memory of the vector flux 
subroutine of the quadrant being calculated. 
It is necessary also to furnish the program at 

1) a9, the value of +(2pz1).0- This value is used in the 


. A 
calculation of _J 


2 
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2) rl, the value of the first trial radius of core 
3) r3, the value of the second trial radius of core 
For a three group solution, values of each of the thirty-six 
fluxes and the three sources at each point are assumed and placed on 


tape in the order indicated in Table I. 


Table I 


Sequence of fluxes at a point 


« i L 
0 0 0 
a 

2 

i 0 

1 

2 

2 0 

1 

2 

1 0 0 
Of < 3 O<i <2 Ox<f <2 
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The points are ordered as indicated in Table II. 


Table II 


Sequence of points 


J m 
0 re) 
i 
2 
pz2 +1 
3 0 
p22 +1 
yy : 
z 
° g. 
@ 
pzl+l1 pz2+l 


As explained in the Procedure the probability of scattering 
from group to group and fission cross section are associated with the 
grid points. These are also placed on the tape of data immediately 
following the flux for f = 3, i = 2, and {= 2 in the order 


indicated on the following page in Table III. 
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Table IIil 


Sequence of additional data for a point 


The total scattering cross section is associated with the region 


between points and are ordered as in Table IV. 


Table IV 


Sequence of Total cross sections 


al m 
0 0 
1 
p22 
2 © 
p2l 0 
° pz2 


Be); 





2 
fe TAPE 110-316-1032 S AND GAUSS CAMPBELL AND PANCIERA 


(2,6) 
p2zle5 pz2=5 pz3=20 = pz y=20 pf1=5000 pf2=10000 pfr3=6)) 
pfhe36 pf5=70), pfl0=2650 pf13=1),70 pf15=300 
/Constants 
al,+1.0 =0.333333 +0.6666667 a2,%0.0 +1.333333 +2.0 
23,4206 aby-1.0 25,41.967705) +1.1306938 
41.306938  +1.967705h a6,+0203 +0.03 +10.0 a7, 40.3183 
40.363); 40.3183 a8,+0,6521),52 40.3)78518 #063h7851,8 
40.6521),52 9,422.0 a10,4+).0 a11,4+1.69))))625 +0, 38):33)1 
/Cross Sections 
bl,#1-0 DITTO 12r/b2, +120 DITTO 12r/m6,+0 m7,+0 
m8, pf3 m9,pf, ml0,+92 mil1,46 ml2,pfl m13,pf2 mih,+0 
m15,+0 
/RADII 


63,4020 rh,+0.0 15,4000 

/NORMALIZING CONSTANTS 

d1}+0.0 d2,40.0 d3,+0.0 rl,40.0 r2,+0.0 810,+0.0 
/NARIABLES 


C1lothe8B  c25"5.6 c3,+0.6885 4002295 © +061590 +3.0355 
+1,0118 +2,0236 

/j-1,m=-1 numbers 4 

nl,#1,0 DITTO 72r/s1, +120 +10 41.0 #1, 41.0 
#100 +120 +10 +1.0 +1.0 +1.0 ' 
/4-1,m numbers 

n2,+1.0 DITTO 72r/s2,41.0 . +1,0 +1.0 t2,41.0 
+1,0 41.0 +1,0 +1.0 41.0 +1.0 
/j-1,m*+1, numbers 

n3,+1.0 DITTO 72r/s3,+1.0 +1.0 +1.0 t3,+1.0 
+1.0 +10 +1.0 +1.0 +1.0 +1.0 
/jy m-1, moumbers 

nl, +1.0 DITTO 72r/sh, +10 +1.0 +1,0 th,+1.0 
41.0 41.0 +1.0 +1.0 +1.0 +1.0 

/jom numbers 

n5,+1.0 DITTO 72r/s5,+1.0 +1.0 +1.0 t5,+1.0 
+1.0 +1.0 +120 +1.0 +1.0 +1.0 

/jom+1, rumbers 

n6,+1.0 DITTO 72r/s6,+1.0 +1.0 +1.0 t6,+1.0 
+20 41.0 +1.0 +1.0 41.0 +1.0 
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/j+lom-1, numbers 


n7,+1.0 DITTO 72r/s7,+1.0 +1.0 +1.0 t7,+1.0 
+1,0 +1.0 +1.,0 +1.0 +1.0 +1.0 

/j+lom, numbers | 

n8,+1.0 DITTO 72r/s8,+1.0 +1,0 +1.0 t8,+1.0 
+1.0 +10 +1.0 +1.0 +1.0 +1.0 

/§+1.m+1, numbers 

n9,+1.0 DITTO 72r/s9, +1.0 +1.0 +1.0 t9,+1.0 
41.0 +10 +10 +1.0 +1,0 +1.0 
/CONVERGED GAMMA 

ml, itaml8 itsr5 idva9 itscel iscO icr3 iscl 
icr6 isce icr2 m2,icac2 idvcl isc2 idvall+c 
iscO imral+c icsl itsc3+c iscO ictm2 isc2 
ictm2 iscO icrpz3 m3,icaal itssl0 icacl itsr1 
* Asel icrpzl+1 isc2 icrpz2+1 ml,ispil isc2 ictml; 

icarl itsr2 iadcl iadcl itsrl iscl ict, 
iscO ictm5S iSTOP m5,icasl10 icpm3 m18,ispo 

/jom SENSE AND CONTROL 

il,itai7 iscl itim6 isc2 itim7 OUT cam7 
mhrLO s1hi5 adm12 tsm15 cam8 mhm6 slhl5 
tsmi) adm15 tsi6 camlL5 suml0O tsm1L5 adm) 
tsm29 caml} sum8 adnk5 tsm26 caml) adm8 
admi5 tsm30 cam7 sul mhm11 slh1l5 atim13 
tsmL5 cam9 mhm6 slhi5 tsm1), adm15 tsm32 
camL) sum9 adm15 tsm31 cam6 dmO cpi2 
cam7 dmO cpi3 IN ispm23 ispm27 ispi5 
i2,cam7 dmO cpil; IN ispm23 ispm21 ispi5 
i3,IN ispm23 ispm2) ispiS i), IN ispm23 ispm16 
i5,iD0B nS 16,+0 +78 i7,ispo 

/DATA READ IN PROGRAMME 

‘m23,itam33 iDIB nl m26,+0 +276, DIB ny 

m29,. +0 276 iDIB n? m30,+0 +276 iDIB 
bl m31,+0 +6 iDIB b2 m32,+0 +6 
m33,ispo 

/ORIGIN PROGRAM 

wL6,itaml7 ispp1l0 ispp13 icaal idvt7 itsdl, ispe3 
iDIB pf13 +2050 +175 ispfl m17,isp0 
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/CENTERLINE PROGRAM 


ispe3 
pfi3 


ispe3 
pf13 


ispe3 
pf13 
+125 
m28 ,ispO 


itstl 
imrt5+2 
imet5+); 
imrt5+10 
Program 
itss5+2 
isusS 
ispe5 


p2,icaa7 
icaa7+); 


imrn5+18 +c 


icaa7+) 
icaa7 
iadt2 
itstl 


iadtl 


p5,icaa8 
icaa8+), 


m21,itam22 isppl, isppl13 
4175 ispfl iDIB 
m22,lspO 

/MIDPLANE PROGRAM 

m2h,itam25 ispp7 ispp13 
+175 ispfl | iDIB 
m25,ispQ 

/GENERAL POINT PROGRAM 
m27,itam28 isppl ispp13 
4175 ispfl iDIB 
iDIB pfl3 +2350 
+2):75 +175 ispf18 
/Source S°,S',S8* at all Points 
pl3,itaplLl icat6 imrt5 
jadt1 itst7 icat6 
fadtl itst7+2 : icaté 
jadtl.  itstl icat6+); 
/Source Normalization+Tolerance 
63 itae7 icat7+2 imrdl 
icat7 imrd] itst7 
icpe5 e6,icsa6 itssl0 
e5,icat7 itss5 e7,ispO 
/SCALAR Flux-General Point 
pl,itap3 isc3 ier3 
imrn5+6%e iactl itstl 
imra8 itstl icaa7 
imrn5+2h+c iadt2 itst2 
imra8 +2 jadti itstl 
ieaa7+2 imrn 5+)2+c 

Ladt2 imra8+), iadtl 
itst2 icaa7+2 imrn5+60+c 
imn5+66+e iadt2 imra6+6 
/Sealar Flux-Centerline Point 
ph, itap6 isc3 ier3 
imrn5+18+c¢ iadtl itstl 
itstl icaa8t+6 = imrn5+5h+c 
icaab imrn5+6+e itst2 


icaa8+2 
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iDIB pea3 
+2h75 +175 
iDIB pf13 
+2225 +125 
iDIB pf13 
+2225 4125 
ispf27  iDIB 
icaa3 imrt5+12 
itstl icat6+2 
itstl icat6+2 
ladtl itst7+) 
icat7+); imrdl 
icpel idvt7 
elidvt7 iadaé 
imrn5+c itstl 
imrn5+12+c 

itst2 .- 
imrn5+30+¢ 
immn 5+36+¢ 
itst2 icaa7+h 
icaa7 imrn5+5)\+e 
iact2 itst2 
itst6+c ictp2 
imrnSt+c itstl 
imrn5+36+c 
ladtl imra7 
imrn 5+2);+c 


+2050 
ispf18 


+2050 
ispf10 


+2050 
ispf10 
pf13 


imrt6+), 
imrt5+6 
imrt5+8 
plh,ispO 


itssS+) 
isua6 
icpe6 


icaa7+2 
Ladtl 
icaa7+2 
iadt2 
itst2 

im 5+1;8+e 


tcaa7t+) 
p35ibp0 


icaa8b+2 
jadtl 
itstl 
iadt2 





itst2 icaa8+l,  imrn5+l2+c iadt2 itst2 icaa8+6 
imrn5+60+c iadt2 imra7+2 iadtl imra2+), itst6+c  ictp5> 
p6,ispO 

/Scalar Flux-Midplane Point 

P7 ,itap9 isc3. icr3 p8,icaa7 imnS+c itstl icaa7+2 
immn5+64+c iadtl itstl icaa7+) imrn5+12+c jadtl 
imra8 itstl icaa7 imrn 5+18¢ itst2 icaa7+2 
imn5+2li+e iLadt2 itst2 icaa7+h imrn5+30¢c iadt2 
imra8+2 iadtl imra2+ itst6+e  ictp8 p9,ispO 

Scalar Flux - Origin 

plO,itapl2 isc3 icr3 pll,icaa? imrn5+c itstl icaa7+2 
immnS+6+e iadtl imra8 itstl icaa7 imrn5+18+c 

itst2 icaa7+2 = immn5+2)i+c jadt2 imra8+2 iadtl 
imra2+) imra2+) itst6+ce ictpll pl2,isp0 

/RADIUS EXTRAPOLATION ROUTINE 

qi ,itag9 icad2 isual itstl icpq2 q3,itst3 icad3 
isual itst2 icpqh q5,itsth isut3 icpq6 icad2 
iexd3 itsd2 icar3 iexr); itsr3 icat3 iexth; 
itst3 icatl iext2 itstl q6,icat3 isua6+2 icpq7 
q8,icad2 isud3 itst6 icar)) isur3 imrtl idvt6 
iadr3 itsr3 q9,ispO q2,imral ispq3 qu,imral ispq5 
/Control 

xloicarl isch icrpz ispml icadl itsd3 x2,icar3 
ispml icadl itsd2 isch ictx3 iSTOP x3,ispql 
ispx2 Q7,iDIB pf5 pf10 pf15 ql0,isppf5 

iSTOP 

/WORKING FLUX 

1470/y1,+0 DITTO 175r/y2,+0 

/First Quadrant Flux 

(Firat Flux 

DA2050/ 1),70/ fl,itaf3 isc3 icr3 f2,icsb2+6+c 
impa5S imrel itstl jadal imn9+c itst2 icas5S+c 
iads9+c imraS imrcl iadt2 imrc3 itste icstl 
inre3 jadal itstl icaal isuc3 imrn6+ce iadt2 
idytl itsn5+c 

/Seeond Flux 

flyicsa2+2 idvrl itstl icsb2+6+c imras itst2 iadtl 
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imrel itst3  iadal+2 imrn9+6+c itsth  icbtl § iadt2 
imecl] itst6 jadal+2 imren9+c iadt) itsth icat6 
isual+2 f5,imrn5+c iadth itst) icasS+tc iads9+c 
imra5 imre] imra2+l, iadth itsth icsal+); iadal+2 
idva2+), itstl ican6+te iadnSte isun9+c isun8+c imrtl 
iadth; imre3+2 itst f6,icst3 imrc3+2 iadal+2 itst2 
icaal isuc3+2 imralt2 itst3 ican6+6+c iadn6+c isund+c 
imrt3 iadth idvte itsn5+6+c 

/Third Flux 

£7 ,icsb2+6+c imraS+2  imrel itstl iedal § imrn9+18+c 
itste icasS+c iads9+c imrad5+2 imrcl iadt2 imre3+6 
itst2 icstl iadal imre3+6 itstl icac3+6 isual 
itst3 ican8+18+c isun9+18+c imrt3 iadt2 
idvtl itsn5+18+c 

/Fourth Flux 

f8,icsa2+2 idvrl itstl icsb2+6tce imra5+2 itst2 iadtl 
imrel itst3 fadal+2  imrn9+2h+c itsth icstl 
iadt2 imrel itst6 jadal+2 imrn9+18+c jadth 
itsth icat6 isual+2 imrn5+18+c jadth itsth 
£9,icas5tc iads9+c imra5+2 imrcl imra2+), iadth itsth 
icsalth fadal+2 idva2+) itstl ican6+18+c iadn5+18+c 
isun9+18+c isun 8+18+c imrtl iadth imre3+8 itstl 
iest3 imrc3+8 iadal+2 itst2 icaal isuc3+8  imral+2 
itst3 ican6+2)+c iadn6+18+e isun5+18+c 

imet3 iadth Ldvt2 itsn5+2lit+c ictf2 £3,isp0 
/Second quadrant flux 

/First flux 

DA2225/ 1170/ £10 oitafl7 isc3 icr3 fll,icsa2+2 

idvr2 itstl icsbl+6+c imra5 itst2 iadtl imrel 
itst3 jadal+l imrn3+12+c itsth iestl jadt2 
imrel itst6 jadai+) imrn3+6+c iadtl itsth icat6 
isual+) imrn5+6+c fl2,iadth itstl iecasStc § iads3tc 8 imrad 
imrcl imra2+) iadth itstl, icaal+) isual+2 idva2+) 
itstl ican6+6+c iadn5+6+c isun3+6+c isun3+6+e itst3 ican6+12+e 
jadn6+6+c iswm5+6¢c imrt3 iadth Ldvt2 itsn5+12+c 
/Second Flux 

fll icsa2+2 idvr2 itstl icsbl+6+c imra5+2 itst2 iadtl 
imrel itst3 jadal+) imrn3+30+c itst) icstl 
iadt2 imrel itst6 fjadal+)  imrn3+2))+c iadth 
itsth icat6 isual+), imrn5+2)+c f15,iadth itsth 
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icas5+c jads3+c imraS5+2 imrcl imra2+); iadth itstl 
feaalt) isual+2 idva2+l itstl ican6+2l;+¢ iadn5+2)+c 
isun3+2h¢c isun2+2)+c inrtl iadth imrec3*+10 itstl 
icst3 jadal+); imrc3+1l0 itst2 £16,icac3+10 isual 
imeal +), itst3 ican2+30+c iadn2+2)+c isun3+30+e 
isun3+2)+c imrt3 Ladtl idvt2 itsn5+30+te ictfll 

f17 ,ispO 

/Third Quadrant Flux 

/First Flux 

DA2350/ 170/ £27, itaf33 isc3 icr3 £28 ,icsa2+2 

idvr2 itstl icsbl+c imra5+l; itst2 iadtl imrel] 
itst3 jadal+); imrn]+)8+c itsth, icstl iadt2 
imrel itstl fadal+h  imrnl+),2+c jiadth itsth 
icatl isual+)  imrn5+l2+c iadth itsth £29 ,icas5te 
iadsitc imraS+); imrcl imra2+), iadtl itsth icaal+) 
isualt2 idva2+l, itstl icanS+h2+c Ladnl+)2+c 
isun2+l2+e isunl+)2+c imrtl iedth imre3+10 itstl 
Lest2 jadal+l) imrc3+10 itst2 icac3+10 isual imral+); 
itst3 ican2+]8+c iadn2+),\2+c isunl+)8+c 
isunl¢+2+c imrt3 ~ iadtl idvt2 £30, itsn5+),8+c 

/Second Fiux 

£31l,icsa2+2 idvr2 itstl icsbltc imra5+6 itst2 iadtl 
imeel itst3 jadal+), imrnl+66+c itstl iecstl 
iadt2 imrc] itstl jadal+),; imrni+60+c iadth 
itsth icatl isual+), imrn5+60+c fadth itsth 

£32 ,icas5+e iadsl+c imra5+6 imrel imra2+), iadth itsth, 
icaal+l, isual+2 idva2+l) itst2 ican5+60+c jadnl+60+c 
isun2+60+e isunl+60+c imrt2 iadth = imre3t+) itstl 
icst3 imrc3+l,; iadal+) itst2 icaal isuc3+), imral+th; 
itst3 icanl+66+c iadnl,+60+c isun5+60+c 

imrt3 iadtl idvt2 itsn5+66+c ictf28 £33,isp0 
/Fourth Quadrant Flux 

/First Flux 

DA2;75/ 170/ £18, itaf26 isc3 icr3 f19,icsb2+c 

impadS+); imrcl itstl iadal inrn7+36+c itst2 
icasSt+e iads7+e imeaSt+); imrcl iadt2 imre3+6 itst2 
icstl iadal imre3+6 itstl icac3+6 isual itst3 
ican8+36+¢ isun7+36+c imrt3 jadt2 idvtl itsn5+36+c 
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/Second Fiux 


f20,icsa2+2 idvrl itstl icsb2¢c imraS+l itst2 iadtl 
imrel itst3 jadal+2 imrn7+h2+c itsth, icstl 
iadt2 imrcl itst6 fadal+2 imrn7+36+c iadth, 
itstl, icat6 isual+2 imrn5+36+c iadth itstl, 
f2l,icasS+c iads7+c imrad5+) imrel imra2+), iadth itsth 
icaalt), idva2+) itstl ican5+36+c iadnl+36+c 
isun8+36+c¢ isun7+36+e imrtl jadth imre3+8 itstl 
icst3 imec3+8 jiadal+2 itst2 f22,icaal isun3+8 imral+2 
itst3 icanhth2+c iadnl+36+c isun5+36+c 

imrt3 jadtl idvt2 itsn5+)2+¢c 

/Third Flux aur 
£23,icsbétc imra5+6 imrel itstl iadal imrnf+5+c 

itst2 icasS¢c iads7+c imrasS+6~ = imrcl iadt2 imre3 
itst2 iestl imre3 jadal itstl icaal isuc3 
imenh+5i+e iadt2 idvtl itsn5+5h+c 

/Fourth Flux 

f2h,icsa2+2 idvrl itstl icsb2+c imra5+6 itst2 iadtl 
imeel itst3 jadal+2 imrn7+60+¢c itsth icstl 
Ladt2 imrcl itstl jadal+2  imrn7+Sh+c jadth 
itsth icatl isual+2  imrnS+5);+c jadt,h § itsth 
f£25,icas5+c iads7+c imra5+6 imrel imra2+¢l, iadth icaalt2 
isual+  idva2+) itstl  icanS+Shtc Vadnh+Sh+tc 
igun8+S)4c isun7+5h+c imrtl jadth imre3+2 8 itstl 
icst3 imee3+2 iadalt2 itst2 icaal isuc3+2 imralt2 
itst3 icanl+60+c Ladnl+5i+e isun5+5l+c 

imt3  iadti idvt2 itsn5+60+e ictf19 £26, isp0 
DA2650/ 70h/ zi, itazlh iscl i erpz1+2 isc2 ierpz2+2 
ieadl isp136213 iscl Z212,itim6 ise2 Z1l,itim7 OUT 
can7 mihi. 0 slh15 adml2 tsmL5 cam8 mhm6 
slhis adm15 tsz2 cam6 dmO cpz3 cam7 
dmO epzh IN ispz5 isppl ispz6 23,cam7 
dmO Cp27 IN isp25 isppl ispz6 Zh,IN 
isp25 ispp7 isp26 27, IN ispz5 ispp10 ispz6 
isc2 ictzll iscl ictz12 Zih,ispO 25,itaz8 iDIB 

nS Z2,+0 +72 z8,ispO z26,ispz6+1 Z1h,ispo 
Zi0,icat6 ispl39z13 icat6é+2 isp139213 icat6é+h ispl36213 29,isp27+h 
213. 
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